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We present an alternative way of tracing the existence of a scalar field based on the analysis of 
the gravitational wave spectrum of a vibrating neutron star. Scalar-tensor theories in strong-field 
. . . gravity can potentially introduce much greater differences in the parameters of a neutron star than 

nJ \ the uncertainties introduced by the various equations of state. The detection of gravitational waves 

, from neutron stars can set constraints on the existence and the strength of scalar fields. We show 

»' ' that the oscillation spectrum is dramatically affected by the presence of a scalar field, and can 

provide unique confirmation of its existence. 
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INTRODUCTION 



(N 

> . 

\^ ' A natural alternative/generalisation to general relativity is the scalar-tensor theory in which gravity is mediated by 
long-range scalar fields in addition to the usual tensor field present in Einstein's theory Scalar-tensor 
theories of gravity can be obtained from the low-energy limit of string theory or/and other gauge theories. The 
existence of scalar fields is crucial in explaining the accelerated expansion phases of the universe e.g. inflation and 
■ quintessence. They are viable theories of gravity for a specific range of the function that couples the scalar field to 
' gravity. Still it is not clear how the scalar fields couple to gravity 0] , a basic assumption is that the scalar and 
gravitational fields Lp and g^^ are coupled to matter via an "effective metric" g^i, = A^{ip)g^^. The existence of 
the scalar field has not yet been verified while a number of experiments in the weak field limit of general relativity 
r ' disproved or set severe limits in the existence and the strength of the scalar fields |^ |3 . 
^jpj' The Fierz-Jordan-Brans-Dicke Q, |^ theory assumes that the "coupling function" has the form A{lp) = aQLp 
• • , i.e. it is characterized by a unique free parameter a§ — (2ti;BD + 3)~^ and all its predictions differ from those of 
. ^ ' general relativity by quantities of order ^3 ■ Solar system experiments set strict limits in the value of the Brans- 
k> " Dicke parameter com i.e. wbd ^ 40000 which suggests a very small < 10~^ see 0, 0, 0, 0. Damour and 
Esposito-Farese [l^, [3 showed that the predictions of scalar-tensor theories in the strong field might be drastically 
[ different from those of general relativity. By studying neutron star models in a simplified version of scalar tensor 
theory where A{ip) = ao<p + /3<p^/2 they found that for certain values of the coupling parameter /3 the stellar models 
develop some strong field effects which induce significant deviations from general relativity. Their claim was based 
on the difference that one can observe in the properties of neutron stars by the introduction of a scalar field (even if 
the coupling constant ao is very small). Damour and Esposito-Farese described the sudden deviation from general 
relativity for specific values of the coupling constants as "spontaneous scalarization" . Harada studied in more 
detail the stability of non-rotating neutron stars in the framework of the scalar-tensor theory and he reported that 
"spontaneous scalarization" is possible for /3 < —4.35. In other words the "spontaneous scalarization" of Damour and 
Esposito-Farese suggests that weak-field experiments cannot constrain the effect of the scalar fields for the strong- 
field regime and prompts for alternative measurements. Such measurements can be based, for example, on accurate 
estimation of the orbital decay of binary systems [3, , in the accurate monitoring of gravitational waveforms from 
neutron stars spiralling into massive black holes |17| or by direct observation of monopolar gravitational radiation 
during the collapse of compact objects. Collapse simulations have shown that indeed a scalar field can be observed by 
LIGO/EGO for the specific range of values of /3 if the event takes place in our Galaxy E [H S El IH IH , while 
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the space gravitational wave detector LISA can also provide constraints in the existence of scalar field 0, |^ . 

Recently, DeDeo and Psaltis |25| suggested that the effects of the scalar fields might be apparent in the observed 
redshifted lines of the X-rays and 7-rays observed by Chandra and XMM-Newton. Testing strong gravity via elec- 
tromagnetic observations is a really novel idea, which have been recently extended, by the same authors p^ . in 
suggesting tests via the quasi-periodic oscillations (QPOs). 

In this paper we examine whether gravitational wave observations of the oscillation s pec tra of neutron stars can 
provide an alternative way of testing scalar-tensor theories. The aforementioned results |lCl lisl show that for a 
specific range of values of the coupling constant (3, which are not constrained by the current experimental limits in the 
weak-field regime, neutron stars can have significantly larger masses and radii. This immediately suggests that the 
natural oscillation frequencies of the neutron stars will be altered accordingly and a possible detection of gravitational 
waves from such oscillations will not only probe the existence of the scalar field but might provide a way of estimating 
its strength. 

The estimation of the stellar parameters (mass, radius and equation of state) via their oscillation properties is not 
a new idea. Helioseismology and asteroseismology are established fields in Astronomy and there is already a wealth 
of information about the interior of our sun and the stars via this approach. In the late '90s it was suggested [23,123 
that the oscillation spectra of neutron stars can reveal in a unique way their properties. The radius, the mass and 
the equation of state can be easily deduced by an analysis of the oscillation spectrum of the /, p and w-modes 29]. 
Moreover, the rotation period of a neutron star can be revealed by the r-mode oscillations since their frequency is 
proportional to the rotation rate (t'.'i— mode ^ 4r2/3). Features such as superfluidity j3Q|| or magnetic fields [sj can 
reveal their presence via detailed analysis of the gravitational and electromagnetic spectra. Recently, it has been 
suggested that compact stars with exotic equations of state, such as strange stars, have a spectrum which carries in 
a unique way their signature [3^ l33j | . 

The oscillations of a neutron star in the scalar-tensor theory will produce not only gravitational but also scalar 
waves m . Detectable scalar waves will be a unique probe for the theory, but this might not be the case if the radiated 
energy is small. Still, we show that there is no need for direct observation of scalar waves, the presence of the scalar 
field will be apparent in the gravitational wave spectrum, since the spectra will be shifted according to the strength 
of the scalar field. 

The paper is organized as follows. In the next section we present the basic equations for the construction of the 
unperturbed spherically symmetric stellar models. We also show the effect of the scalar field on the stellar structure. 
In Section 131 we derive the perturbation equations which will be used for the numerical estimation of the oscillation 
frequencies. Finally, in Section 0] we discuss the results in connection to gravitational wave asteroseismology. Note 
that the present study will be based on two equations of state, which have been used earlier [lOllTsl . Our aim is to 
demonstrate the effect of the scalar field in the neutron star oscillation spectra, while more detailed analysis for a 
wide range of EOS is underway. 

2. STELLAR MODELS IN SCALAR-TENSOR THEORIES OF GRAVITY 

In this section we will study neutron star models in scalar-tensor theory of gravity with one scalar field. This is 
a natural extensions of Einstein's theory, in which gravity is mediated not only by a second rank tensor (the metric 
tensor g^^^) but also by a massless long-range scalar field (p. The action is given by 0] 

S = J V^* {R* - 2.9r^,p¥',.) d^x + 5™ [^m, A^{ip)g,^,] , (1) 

where all quantities with asterisks are related to the "Einstein metric" g^,fi,^, then i?* is the curvature scalar for 
this metric and G, is the bare gravitational coupling constant, represents collectively all matter fields, and 

Sm denotes the action of the matter represented by ^',„, which is coupled to the "Jordan-Fierz metric tensor" g^i^. 
The field equations formulated better in the "Einstein metric" but all non-gravitational experiments measure the 
"Jordan-Fierz" or "physical metric" . The "Jordan-Fierz metric" is related to the "Einstein metric" via the conformal 
transformation, 

Ot^iy = A^{'^)9*tJ^u- (2) 

Hereafter, we relate all tilded quantities with "physical frame" and those with asterisk with the "Einstein frame". 
From the variation of the action S we get the field equations in the Einstein frame 
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where T!^'^ is the energy-momentum tensor in the Einstein frame which is related to the physical energy-momentum 
tensor T^^ as follows, 



The scalar quantities and a{ip) are defined as 



SS„ 



9^ 



dlnA(^) 



a{Lp) = 



dip 



(5) 

(6) 
(7) 



It is apparent that a{ip) is the only field-dependent function which couples the scalar field with matter, for a{^p) — 
the theory reduces to general relativity. 

Finally, the law for energy-momentum conservation \/ ,jT^^ = is transformed into 



(8) 



and we set (^q the cosmological value of the scalar field at infinity. In this paper, for simplicity, we adopt the same 
form of conformal factor A{ip) as in Damour and Esposito-Faresejlflj . which is 



(9) 



i.e. OL(ip) = jiip where /3 is a real number. In the case for /3 = 0, the scalar-tensor theory reduces to general relativity, 
while the "spontaneous scalarization" occurs for /3 < —4.35 [l5l |. 

We will model the neutron stars as self-gravitating perfect fluids, made out of degenerate matter at equilibrium and 
admitting cold equation of state. Then the metric describing an unperturbed, non-rotating, spherically symmetric 
neutron star can be written as 



dsi = g^^^dx^'dx" = -e^'^dt^ + e^'^dr^ + r\de^ + sin 



where 



-2A 



Mr) 



(10) 



(11) 



while the "potential" function (f>(r) will be calculated later. The stellar matter is assumed to be a perfect fluid 

T^,u ={p + P] U^U, + P~g^,. (12) 



where C/^ is the four-velocity of the fluid, p is the total energy density in the fluid frame, and p is the pressure. 

Spherical symmetry simplifies significantly the procedure of constructing stellar models. Using equations |(2Jl, Q 
and , we can obtain the following set of equations for the background configuration 0, 0| 



— = 47rG. 



dr 
d$ 
dr 
dip 
dr 
d^ 
dr 
dP 
dr 



^TiG^r^A'p + -r{r - 2/i)*^ 



''A^P 1 



r- 2/1 



r{r — 2/i) ■ 



— = 47rG» 



rA* 



r- 2/i 

p + P 



a(/5 - 3P) + r(p - F)^' 



2(r-M) 
r(r — 2/i) 



dr 



47rG, 



HA^P 1 



r - 2/t 



r(r — 2p) 



(13) 
(14) 
(15) 
(16) 

(17) 
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Near the center the background quantities P and fi can be expanded as, 



ip{r) 


— in -\- 


-I- (^(^r J, 




/ 

= = 


(y52r + 0(r^), 


$(r) 






P{r) 


= Pc + 


ip2r2+0(r4), 




= 0{r^ 





where the expansion coefficients are given by 



2 =47rG,A*P„ 
<P2 = — G,A^ ttc (^/?: 



3P, 



P2 = -[Pc + Pc]{^2+acip2) 



(18) 
(19) 
(20) 

(21) 
(22) 

(23) 
(24) 
(25) 



where Ac = A{ipc) and ac = a{Lpc)- 

Outside the star, the metric is assumed to be static and spherically symmetric and since we deal with only one 
scalar field, it gets the form P| 



dsi = ^e^dt' 



{de^ 



sm 



b/a 



r{r — a), 



where r is a radial coordinate, given by the following relation 



l-b/o 



(26) 
(27) 

(28) 



Moreover, in the above metric (|27|l . a and b are some constants, which are connected with the total scalar charge uja 
and the total ADM mass Mpj^yi ,i.e., — b"^ = 4:Lua^ and b — 2Mjsj^y[. Finally, the asymptotic form, at spatial 
infinity, for the metric and scalar field will be given as functions of the total ADM mass and the total scalar charge 



UJA 



tp = ipo-\ h O 



where ry^jy is the Minkowskian metric. 

By matching this exterior solution and the interior metric, one gets the following relations ^ 



(29) 
(30) 



1 - 



2fis 
R 



1/2 



' 2 ' R{R-2ns)' 



cxp 



zarctanh 



larctanh 



*2 



(31) 
(32) 
(33) 



where the functions with the subscript s correspond to their actual values at the stellar surface and the prime denotes 
the derivative with respect to r. A more general scheme for fixing the asy mptotic value of wq can be derived using 
the cosmological model of Damour and Nordtvedtjs^l, see for example |33- Still, this approach has been applied for 
/3 > which is not the subject of our present study. 
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In order to determine the stellar properties, an additional equation is needed, i.e. the equation of state (EOS). Here 
we will use a polytropic one given by 

~ TL 

P ^ Knomb — , (34) 
no 

P 

p^nmb + Y^—j, (35) 

rrib = 1.66 x 10^24 (3g) 

no = 0.1 fm"^ (37) 

We have selected F = 2.46 and K = 0.00936 in agreement with a fitting to tabulated data for EOS A 34], and 
r = 2.34 and K = 0.0195 to fit EOS II 35]. In other words the parameters T and K are adjusted to fit the tabulated 
data from these two realistic equations of state which have been used in earlier studies of the problem ^0, 15]. The 
present study will be constrained to these two EOS since our aim is to demonstrate the effect of the scalar field in 
the neutron star oscillation spectra while more detailed analysis for a wide range of EOS is underway. 



2.1. Neutron Star Models 



Here we present some typical stellar models for the two equations of state under discussion i.e. EOS A and EOS 
II. For the construction of the stellar models, there exist three freely specifiable parameters, these are: the constant /3 
of the conformal factor, the value ipo of the scalar field ip at infinity, and the central density pc- We will also consider 
positive values of ipo, because the basic equations (|13|) - H17() are symmetric under the reflection (p —> —(p. Also, we 
will only consider stellar models with f3 = —6 and /? = —8, for these range of values the effect of the scalar field is 
more pronounced. Binary pulsar data already suggest larger values for /? i.e. f3 > —4.5, but this work is based on a 
different framework and the results can be used as an alternative way of constraining the appearence of "sponatancous 
scalarization" in neutron stars. 

In Figuren]we show stellar models with (3 = —6. In the left column we plot models for the EOS A and in the right 
column models for the EOS II. In the two upper panels we plot the ADM mass Af^^ji^ as functions of the central 
density, in the two middle panels we plot the value of the scalar field in the center of the star, ipc as functions of 
the central density. It is apparent that the effect of the scalar field is pronounced for central densities higher than 
5 X 10^'*gr/cm'^, for lower densities there is no way that one can trace the existence of a scalar field. The lower panel is 
a "typical" mass-radius relation where one can observe the dramatic effect of the scalar field in the stellar structure. 
For comparison in every panel the general relativistic stellar models are shown by a solid line (GR) 

In Figure 121 we present stellar models for f3 = —8, the ordering of the graphs is as in Figure ^ Now the effect of 
the scalar field is more pronounced than before. Its presence can be traced for even lower densities while for higher 
densities the stellar models divert considerable from those of GR. The effects on the mass and radius are even more 
dramatic in the mass-radius diagrams in the lower panels. 

Finally, in order to stress the effect of the scalar field on the maximum mass of neutron stars we show in two tables 
the relevant parameters for each EOS. The maximum masses are well beyond the observed values of neutron star. 
This is also true for the observed redshifts. It is worth mentioning that the redshift of the maximum mass models is 
quite unaffected by the changes in the central density and the scalar field parameter. Although for these maximum 
mass models the redshift is unusually large, for the typical models that we will use the surface redshift z is below the 
maximum observed redshift which is 2; = 0.35 38]. This means that neutron stars of the scalar-tensor theories cannot 
be exclude by the present electromagnetic observations. 



TABLE I; The stellar parameters for models with maximum ADM mass (EOS A). 



/3 


tpo 




pc (g/cm^) 


R (km) 




Mat)m/R 


z 


0.00 




1.655 


3.742 X 10^^ 


8.56 




0.285 


0.527 


-6.00 


0.00 


2.054 


3.048 X lO^'^ 


11.64 


0.2676 


0.261 


0.445 


-6.00 


0.03 


2.462 


3.199 X 10^^ 


13.67 


0.3269 


0.266 


0.462 


-8.00 


0.00 


3.010 


3.294 X lO^'^ 


16.49 


0.3338 


0.2695 


0.473 


-8.00 


0.03 


3.792 


3.317 X 10^^ 


20.61 


0.3755 


0.2717 


0.480 
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TABLE II: The stellar parameters for models with maximum ADM mass (EOS II). 



/3 




a^admA^o 


Pc (g/cm^) 


R (km) 






z 


0.00 




1.946 


2.543 X lO^'^ 


10.90 




0.264 


0.454 


-6.00 


0.00 


2.423 


2.282 X lO^'^ 


14.00 


0.2724 


0.256 


0.430 


-6.00 


0.03 


2.896 


2.423 X lO^'' 


16.32 


0.3307 


0.262 


0.450 


-8.00 


0.00 


3.529 


2.523 X 10^^ 


19.54 


0.3366 


0.267 


0.464 


-8.00 


0.03 


4.442 


2.543 X 10^^ 


24.42 


0.3792 


0.269 


0.470 



3. PERTURBATION EQUATIONS IN THE COWLING APPROXIMATION 

In this section we present the perturbation equations for non-radial oscillations of spherically symmetric neutron 
stars in scalar-tensor theory. We derive the equations, describe the numerical method for getting the spectra and 
finally we present results for a few characteristic values of the scalar field j3. The presentation in this section is quite 
extended, we have chosen to show all the details of the calculation since this is the first article dealing with the stellar 
perturbations in scalar-tensor theory and there is a wealth of new functions and notations. We will also keep arbitrary 
the functional form of the scalar function a{ip) which will be fixed only later during the numerical calculations. In 
this first investigation we will not consider the full problem but instead we will restrict it to the so called "Cowling 
approximation" . This means that the fluid is perturbed on a fixed background, in this way we freeze the spacetime 
and scalar field perturbations i.e. (5.9^,^ = and 5lp = 0. The Cowling approximation limits our study to the modes 
which are directly related to the fluid perturbations i.e the /, p and fl-modes, while we cannot study the emission of 
gravitational and scalar waves neither the families of the spacetime 39] and scalar modes. 

The fluid perturbations will be described by the Lagrangian displacement vector 

where W and V are functions with of t and r. The unperturbed 4-velocity tj'^^^f^ is 

hereafter an index "{By^ will denote the unperturbed quantities. The perturbed 4-velocity (5f/^ has the form 

Jf/M= o,e A _ ^Yi.^^ (40) 



dt ' dt de' sin^ e dt d(j)J A{(p)r^ 
The perturbed energy-momentum tensor 5T^'^ in the Cowling approximation gets the form 

+ (p(^) + P^^yjglfiu^^^'^su'' + 5tj^u^''>) + sPs;;. (41) 

The Lagrangian variation of the baryon number density An is 



Vrr-:TI7k' (42) 



where V). denotes the covariant derivative in a 3-dimensional with metric g^iu and A the Lagrangian variation. 
In this equation, the first term of the right hand side expresses the 3-dimensional divergence of the fluid, and the 
second term is the amount of the volume change due to the metric perturbation, but because we use the Cowling 
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FIG. 1: Stellar models with fi = —6 are shown. The left column corresponds to models for the EOS A while in the right to 
models for EOS II. In the upper panels the ADM mass Mp^y^ and the value of the scalar field in the center of the star, are 
plotted as functions of the central density. The effect of the scalar field is apparent for pc > 5 x lO^'^gr/cm'^. The lower panel 
is a mass-radius diagram where the dramatic effect of the scalar field can be observed. In every panel, the general relativistic 
stellar models are shown by a solid line (GR) 



approximation we neglect this second term. By employing the above perturbation variables, equation H42|) for the 
Lagrangian variation of the baryon number density will be written as 



An 

h 



m 1(1 ^\)v 

dr 



-A 



W 



Yin 



(43) 



For adiabatic perturbations, using the first law of thermodynamics, we can get the following relation between the 
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FIG. 2: Stellar models with /3 = —8 are shown. The left column corresponds to models for the EOS A while in the right to 

models for EOS II. In the upper panels the ADM mass A^ADM ^'^'^ value of the scalar field in the center of the star, tpc are 
plotted as functions of the central density. The effect of the scalar field is apparent for pc > 3 x 10^''gr/cm"^. The lower panel 
is a mass-radius diagram where the dramatic effect of the scalar field can be observed. In every panel, the general relativistic 
stellar models axe shown by a solid line (GR) 



change of the baryon number density and the pressure 



+ P 

Ap = !-—Ah. (44) 

n 



Therefore we can express the Eulerian density variation i5p as, 



5p=(p(-)+p(-))^-^e-^y,„, (45) 
n or 
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here we have used the relation between the Lagrangian perturbation A and Eulerian perturbation S 



(46) 



From the definition of the adiabatic constant 
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ainP\ _ fiAP 
9 Inn / 



PAn 



(47) 



we can derive the form of the Eulerian variation of the pressure 



n 



An dP^^'> 



dr 



rYir 



Finally, the combination of equations H44(l and H47|l. provides the standard form of the adiabatic constant 



(48) 



(49) 



By taking a variation of the equation for the conservation of energy-momentum '^vT^ = 0, we can get the equation 
describing the perturbations of the fluid in the Cowling approximation 



The analytic form of the above equation for the values of the index /i = 1 , 2 is : 
Vjf^" = 



1 



^ {~P + P)—e 



d_ 

dr 



~ , 1 _^dW l{l + l)V ^ ^ _,^W 



dP 
dr 
dPdP 
dr dp 



dr 



1 _[^dW l{l + \)V ^ ^ _f^W 
-776 '^-^ h TT-^ h Safe — 



dr 

}_^-KdW_ ^ l{l + l)V 



dr 



dPW 



-AdpW 



d<S> 
— + a* 
dr 



dr 



Vjf/ = 

^ ip + P)e 



-2$ 



df^ 



1 ^^dW , l{l + l)V , ^ 
— e — 1 2 ^ 3aWe — 



dr r 



dr 



(50) 



(51) 



(52) 



By assuming a harmonic dependence on time the perturbation functions will be written as W{t, r) = W{r)e^'^^ and 
V = V{t, r) ~ V{r)e^'^* and the above system of equations gets the form 



u;\p + P)e^-^'^^ + ± 



-- , 1 .dW lil + l)V ^ ^ .w 



dr 



„-A 



dr 



dPW 

dr r^ 



dP 
dr 

dp dP 
dr dp 



1 .dW l{l + l)V ^ ^ jyW 



dr 



1 _^dW l{l + \)V ^ ^^W 
re — 1 h oaWe —tt 



— + a* e"^^ — 



dr 



j,dpW 



dr 



-u;^e-^''{p + P)V + 



jP 



dr 



dr 



dPW 



dr ) dr r^ 



+ e 



^dPW 
dr 



= 0. 



(53) 
(54) 
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A further simplification, can be achieved, for the first of the above equations by using an appropriate combination of 
the form (j (|54|) /dr — (|53|l . this leads to the equation 



dr 



ip + P)V -uj\p + P)e 



A-2* 



W 



dP 

dr 

dp dP 
dr dp 



1 _j^dW l{l + l)V ^ ^ 



dr 

1 _^dW lil + l)V 



dr 



T 1 -A^pW 



dr 



dr 



+ 



W 



dr 



a* e 



^dPW 
dr 



which can be further simplified by proper substitutions of dP/dr and dW/dr from equations (|17|) and H54|l 



dV d$ 



= 2—V-e" 



W 



(55) 



(56) 



Thus from equations (|54() and H56() . we get the following simple system of couple ODEs 





dp 


dr 


~ dP 








= 2— 




dr 



dr 



/(/ + l)e'^V - Sa'^W, 



W 



(57) 
(58) 



which together with the appropriate boundary conditions at the center and at the surface constitute an eigenvalue 
problem for the parameter lu (the eigenfrequency) . 

The above system of perturbation equations and near the stellar center gets the following simple form 



dW 
dr 

dV W 



dr 



+ ;(? + i)F«o, 

0. 



with the following set of approximate normal solutions 



W{r) = Br^+^ + 



Vir) 



B 

r 



(59) 
(60) 

(61) 
(62) 



where B is an arbitrary constant. The two approximate solutions near the center, Eq. 1)61(1 and ((62(1 . suggest the 
introduction of two new perturbation functions, W{r) = W^(r)r'+^ and V{r) = V{r)rK After this change, the 
approximate solutions near the center become 



Wir) ^W, + -W2r^ + • • • , 
V{r) = K + ^V2r^ + • • • , 
where Wc = B, Vc = —B/l and the coefficients W2 and V2 of the second order terms are 



V2 = 



1 



21 + 3 

W2 = 4$2V; -{1 + 2)^2- 



3aaip2Wc + 2{l + 3)^2Vc ~ ^ {io^e'^'^^Va + ($2 + acip2)Wc} 

P2 



(63) 
(64) 

(65) 
(66) 



At the surface, the boundary condition is the vanishing of the Lagrangian perturbation of the pressure AP = 0. 
The Lagrangian perturbation of the pressure is described by AP = -fPAn/h, and thus making use of equation 1(47(1 . 
we get the following boundary condition at r = R 



ip + P) 



2 -2*T/ , -A 

LU e V + e 



d$ 
dr 



W 



0. 



(67) 
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For this new functions the system of perturbation equations becomes 

1(1 + 1] 



dW 
dr 
dV 
dr 



dp 

Hp 



rV 



d<P I 



dr 



V -e" — 
r I r 



d<P 
dr 
W 



l + l 



+ 3a* W 



with following boundary conditions at the center and the surface 

W = -IV (at r = 0), 
iu^e-^'^rV + e-^(^+a-^\w^O (at r - i?). 



(68) 
(69) 

(70) 
(71) 



Equations H68|l . (|69l) . H7U|I and (|71|) form a well posed eigenvalue problem for the real eigenfrequency w^. Using 
shooting method we can get in a quite simple way the eigenvalues i.e. the characteristic frequencies of the oscillations. 
It should notice that the above system of perturbation equations and boundary conditions has been derived for an 
arbitrary form of the scalar function a{ip), later in the numerical calculations we will fix its functional form to the 
one of equation a{ip). 



16. 
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— 1 1 r- 
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-/3 

FIG. 3: The normalized eigenvalues lu for the first few modes (/, pi, p2, pa and p4) are shown as functions the parameter /3 
for the equations of state EOS A (solid line) and EOS II (dashed line). The eflect of the "spontaneous scalarization" is more 
pronounced for the higher modes. The asymptotic values of the scalar field and of ADM mass are fixed to the values tpo = 0.0 
and Madm = IAMq. 



3.1. The Oscillation Spectra 

The spectrum of an oscillating neutron star is directly related to its parameters, mass, radius and EOS As we 
have seen in Section El the presence of the scalar field on the background star is influencing both the mass and the 
radius, while in the way that it enters in the equilibrium equations it has a role of an extra pressure term i.e. it seems 
to alter the actual EOS. As before we will restrict our study only two equations of state i.e. EOS A and EOS II, and 
we will use models from those described in the previous section [2. II 

First, we will examine the effect of the scalar factor (3 on the frequency and especially whether the "spontaneous 
scalarization" can be traced in the spectrum. A discontinuous change in a system, as one varies its parameters, signals 



12 



a catastrophic behavior and the so called "spontaneous scalarization" [loj is related to it. Harada shown that 
in the scalar-tensor theory the value (i — —4.35 is the critical one for the "spontaneous scalarization". DeDeo and 
Psaltis [23 verified that the effects of the scalar field on the line redshifts of neutron stars become pronounced for 
values fi < —4.35. We show here that the spontaneous scalarization is also present in the oscillation spectra of neutron 
stars. In order to study this effect we fixed the asymptotic value of the scalar field to (^0 — 0-0 and we constructed 
a sequence of stellar models with M = I.4M0 by varying /3. In Figure |2I the effect of varying /3 on the frequencies 
of the fluid becomes immediately apparent. For values (3 < —4.35 the frequency gets a sharp change and increases 
linearly with decreasing (3 signaling the "spontaneous scalarization" . The effect becomes more pronounced for the 
higher p-modes as du}/d{—f3) increases with the order of the mode showing the dramatic impact of the "spontaneous 
scalarization" on the spectrum. By studying the first thirty modes for each equation of state we have actually found 
that: 

" for /3<-4.35 (72) 



a(-/3) 4 

where uJn is the frequency of the nth mode. This relation seems to be independent of the equation of state and 
suggests that a possible tracing (via electromagnetic or gravitational observations) of the higher p-modes will signal 
the existence of a scalar field even for small deviations from general relativity. 

The idea of a gravitational wave asteroseismology was based on empirical relations that can be drawn for 

the relation of the stellar parameters to the eigenfrequencies of an oscillating neutron star. These empirical relations 
were derived by taking into account data for a dozen or more EOS, and it was shown that through these relations 
one can extract the stellar parameters by analyzing the gravitational wave signal of an oscillating neutron star. The 
easiest relation to be understood on intuitive physical grounds is the one between the fundamental oscillation mode 
the /-mode and the average density. This relation emerges naturally by combining the time that a perturbation 
needs to propagate across the star and the sound speed, this boils down to a linear relation between the period of 
oscillation and the average density or the star or better ~ M/B?. This was the reason that we have normalized 
the frequencies in Figure 13 with the average density. According to [2^ the empirical relation between the /-mode 
frequency and the average density of typical neutron stars is 

/._..(kHz). 0.78+ 1.63 (im^) (^3) 

Almost, all EOS follow this empirical relation (including the two EOS used in this article). This observation suggests a 



= 0,0 (GR) 




EOS A 



J 1 I 1 I I I 1 1 1 I 1 I 1 I L. 




EOS II 



J 1 I 1 I L. 



0. 020 0. 026 0. 030 0. 036 0. 040 0. 046 



0. 020 0. 026 0. 030 0. 036 0. 040 0. 046 



, / ^ADM , / ^ADM 

V i?3 V i?3 

FIG. 4: The frequency of the /-mode as function of the averaged density (Madm/^?^)^''^ of the star (notice that // = ujf/2%). 
The thick solid line corresponds to the values of the mode for /3 = 0.0 (GR) while we also show the effect of the scalar field for 
three values of the scalar parameter {3 (-8.0, -7.0 and -6.0). Here we have assumed that <^o = 0.0. The left panel corresponds 
to EOS A and the right panel to EOS II. 
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unique way in estimating the average density of a star via its /-mode frequency, and it can be a very good observational 
test for the neutron stars in scalar-tensor theories. In figure 0] we draw the frequency as function of the averaged 
density. The actual relations of the frequency of the /-mode as functions of the average density for general relativistic 
neutron stars are shown as thick solid lines in both panels of figure 01 Both solid lines correspond to normal neutron 
stars with f3 = and follow the empirical relation (|73l) with quite small error. The introduction of a scalar field 
even in moderate central densities alters completely the behavior of the /-mode for both EOS. The frequencies grow 
considerably faster as functions of the average density, for /3 < —4.35, and the three examples that we have chosen 
(/3 = —6.0, -7.0 and -8.0) show exactly this behavior. The change is quite dramatic even for typical neutron stars 
with average density. Depending on the value of the parameter /3 they become 30 — 50% larger than those of a general 
relativistic neutron star. This can be an observable effect, since the detection of frequencies are higher than those 
expected from a typical neutron star, would signal in a unique way the presence of a scalar field. 




4. 



1.0 2.0 3.0 4.0 1.0 2.0 3.0 4.0 

FIG. 5: The normalized frequencies of the first four fluid modes / (upper-left panel), pi (upper-right panel, p2 (lower-left 
panel), and ps (lower right panel) are plotted as functions of the ADM mass. The results of the two EOS are shown for three 
values of the scalar parameter (3 (0,-6 and -8). The asymptotic value of the scalar field is ipo = 0. 

A careful study of Figure El suggests that the observed higher frequencies predicted by scalar-tensor gravity might 
be attributed to a denser neutron star. This possibility cannot be ruled out if only one mode is observed, but as one 
can easily find out from Figures 01 and [3 the effect of the scalar field on the higher fluid modes, e.g. pi, p2 etc are 
more dramatic, and a synchronous observation of a few modes will not only probe the presence of a scalar field but 
it might provide a direct estimation of its strength. Concluding the discussion related to Figure 01 we should admit 
that there is no apparent explanation why the /-mode frequencies after reaching a maximum either do not increase 
as the density increases (/? = —6) or even move towards lower frequencies (/? = —7.0 and -8.0). 
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In the two preceding Figures 13 and 0] we have shown the effect of the varying scalar factor /? on the frequency of 
the modes. In Figure O the ADM mass and the asymptotic value of (po were fixed while in figure 01 we have studied 
only one mode (the /-mode) for fixed ipQ. The next figures show the effect of the scalar field on the frequencies of 
the /, pi, p2, and modes for varying ADM mass or/and the asymptotic value of ipa. In Figure [S] we plot the value 
of the normalized frequency of the mode as function of the ADM mass. One can easily observe that the differences 
due to the presence of the scalar field are really impressive. The models have been chosen to span a range of masses 
from IMq up to the maximum allowed mass from each theory. The normalized frequencies in these diagrams show 
that the imprint of the scalar field is not only apparent but it can influence the spectrum in a dramatic way. For 
example, for all four oscillation modes (/, pi, p2, and p^) the normalized frequency of the maximum mass model in 
GR is half of the corresponding frequency to an equal mass model with /? = — 8 or about 40% smaller for (3 — —7. 
A very interesting observation is that in contrast to the redshift of the atomic lines |25l| which for the maximum mass 
models is typically smaller than 10% (see Tables I and II) the difference in the frequencies between the models with 
/3 = 0(GR) and /3 = —8 is larger than 50%. This figure suggests that a possible observation of more than one mode 
can clearly signal the existence of the scalar field. 

In all previous discussion about the effect of the scalar field on the eigenmode frequencies of the star we have shown 
results for stellar models for which the asymptotic value of the scalar field was assumed zero. It is still an important 
question whether the asymptotic value of if can be traced via the neutron star asteroseismology. The results of a 
varying ipo are shown in Figures IHl [3 Figure IHI is similar to Figure 01 but here we have fixed the scalar factor to 
/3 = —6 and we vary ipQ. In figureHwe show the dependence of /, pi and p2 modes for varying ipQ for the EOS A, the 
results for EOS II are similar. The last two figures show that the imprint of an asymptotically non-vanishing scalar 
field can be observed, while its becomes more pronounced for the higher modes. Varying tpQ from to 0.03 we observe 
frequency variations larger than 10-20%, variations of this order can provide additional constraints on the asymptotic 
value of the scalar field. 




0. 020 0. 025 0. 030 0. 035 0. 040 0. 045 0. 020 0. 025 0. 030 0. 035 0. 040 0. 045 

V i?3 V i?3 

FIG. 6: The frequency of the /-mode as function of the averaged density {Maum/ R^)^^^ of the star (notice that // = aj//27r). 
The thick solid line corresponds to the values of the mode for l3 — 0.0 (GR), the thinner continuous and dashed lines show 
the effect of the asymptotic value of the scalar field tpo for a fixed value of the scalar parameter j3 = —6.0. The left panel 
corresponds to EOS A and the right panel to EOS II. 



4. CONCLUSIONS 

In this paper we have discussed the effect of the scalar-tensor theories on the oscillation spectra of neutron stars. 
Scalar-tensor theories of gravity are generalizations of general relativity and provide a natural connection to superstring 
theories ,^,40j as well as to inflation 41]. The presence of a scalar field in the the strong-field regime has already been 
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FIG. 7: The dependence of eigenfrequencies on the value of ipo, where it is fixed /3 — —6 for EOS A. The range of ADM mass 
is from I.OMq to maximum ADM mass for each ipo- The open circles, squares and triangles represent the /, pi and p2 modes, 
respectively. 



constrained by weak field experiments and here we provide an additional way of testing/constraining their existence 
via the gravitational wave asteroseismology. 

We show that the oscillation frequencies of neutron stars carry very "clean" imprints of the presence of the scalar 
field. An observation of the neutron star oscillation spectrum via gravitational waves or via electromagnetic signals 
emanating from or around the surface of a neutron star will not only probe the existence of the scalar field but it 
might also provide a measurement of its asymptotic value. 

In this study we have used the Cowling approximation which, although a restricted way of studying stellar oscillation, 
has proven to be a very accurate and useful tool in asteroseismology. Still, more detailed study is needed for proper 
modeling of the effect. The inclusion of metric and scalar field perturbations will result to additional information in 
the gravitational wave spectrum. The combination of this extra information should provide more accurate constraints 
on the existence of the scalar field. 
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